Kondo effect versus Jahn- Teller effect: A new class of quantum impurity problems 
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Impurity ions with orbitally degenerate ground electronic configuration embedded in metals are 
subjects to both the Kondo effect and the Jahn- Teller instability. While increasing the Jahn- Teller 
coupling, one passes, from the region of parameters where the perturbation theory is applicable, 
through a complicated intermediate region, to a region where phonons, localized electrons and 
conduction electrons form a coherent state. 



In the course of time, the Kondo model, originally pro- 
posed in an attempt to account for a resistance minimum 
in magnetic alloys, became a kind of reference model 
in condensed matter theory, revealing nontrivial infrared 
properties of metals. It is the simplest, single channel, 
version of the Kondo model that has required decades 
for its theoretical understanding serving as a laboratory 
for applying different methods Jj]]. Along with this de- 
velopment, though somewhat delaying, efforts to devise 
realistic models describing magnetic impurities in metals 
were taken. These investigations date from the discus- 
sion of an impurity ion <i-level degeneracy in the original 
Anderson paper {§ and (after § and a group theoretical 
analysis of |Q) they led to the paper || where a class 
of models, now known as multi channel Kondo models, 
has been sorted out and an attention to a peculiar ( "non 
Fermi liquid" ) low energy behaviour of these models was 
drawn. Bethe ansatz solutions for multi channel Kondo 
models were shortly reported M. During last few years, 
numerous studies of the unconventional low energy fea- 
tures have been carried out. They include large k (k - 
number of channels) expansion conformal field the- 
ory H and bosonization |9) approaches. Despite this es- 
sential success in studying the models, their relevance 
to concrete impurity problems has not been explored in 
detail, basically remaining on the level of Ref. §. A 
remarkable exception to this is given by comprehensive 
investigations of two-level systems in metals whose low 
energy behavior was conjectured to be described by the 
two channel Kondo model JhJ; the conjecture being re- 
cently confirmed by point contact spectroscopy experi- 
ments |pd| . 

One of the main sources of motivations for studying 
multi-channel Kondo models is provided by unusual tem- 
perature (magnetic field) dependence of various physical 
quantities in heavy fermion compounds |]l2| . This issue 
has been raised by Cox for the case of Uranium com- 
pounds fl3|| and, recently, for some Cerium compounds 
fl4j| . At least for dilute systems, when coherence effects 
can be neglected, single impurity properties are likely 
to be dominants (though, even for dense compounds, 
there are considerable empirical evidences favouring local 



physics). Fortunately, new (and puzzling) experimental 
data on dilute U impurities in MPd^ (with M being Y, 
Sc, La, etc.) became recently available fTEfl . Fermi level 
tunning jl6j and ionic radii p"5| arguments evidence U 
being in U* + ionic state. The J = 4 Hund's rule ground 
state is split by the cubic crystal field into Ti, T^, T4, 
and r$ states. Admitting that the point charge model of 
Ref. fl^l is oversimplified, one may assume 

U the T 3 

nonmagnetic doublet being the ground state (in agree- 
ment with neutron scattering data Jl1| though the latter 
experiments are performed on rather heavily doped sam- 
ples). That would lead to a two channel Kondo model 
explaining part of the data - the logarithmic low 
temperature specific heat and the magnetic susceptibility 
behaviour |]l9| but not the linearly (in temperature) sat- 
urating resistivity which remains to be reconciled. This 
line of arguments concerns purely electronic properties 
disregarding lattice degrees of freedom. An important 
point is, however, that an orbitally degenerate ground 
state of the impurity ion realizes classical conditions for 
the Jahn- Teller instability poj . Crucial difference with 
the standard situation is that local phonon modes inter- 
act here not only with localized electrons but also with 
those from conduction band (coupled, in turn, to the de- 
generate local electronic configuration) . Apart from some 
brief comments (e.g., in Ref. p^]), this interesting prob- 
lem has been overlooked in the literature (and poorly 
understood). 

This letter is intended to give a qualitative physical dis- 
cussion of that complicated interplay of the Kondo- and 
the Jahn- Teller effects realized by orbitally degenerate 
ions in metals where the three ingredients of the system 
(the local electronic configuration, the conduction elec- 
trons, and the local phonon modes) are tied together. 
The main trends are explored whereas detailed calcula- 
tions (anyway not allowed by the format of the letter) 
will be presented elsewhere 

The problem exists for any orbitally degenerate ionic 
state (T3, T4, and Tg for the cubic point group), mag- 
netic or not. Low energy properties are of interest, so 
the temperature is understood to be small as compared 
to the crystal field splitting energy (on the scale of 100K 
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for U — MPd.3 if estimated from the specific heat Schot- 
tky anomaly [pUf), Generally, the Hamiltonian of the 
problem can be written in the form: 



H = H P 



ph 



Hel-ph + Hk + HjT, 



(i) 



where H e i describes conduction electrons, H p h accounts 
for the lattice, H e i_ p h represents conventional electron- 
phonon coupling, and Hk and Hjt stand respectively for 
the Kondo- and the Jahn- Teller couplings to be specified 
shortly. 

As an example the cubic point group and the T 3 state 
is considered in this letter. The r 3 representation is 
spanned by localized electronic wave functions transform- 
ing as |0) oc 3z 2 — r 2 and |e) cx x 2 — y 2 . [T§] = ri+r 3 and 
the operators (again spanning r 3 ), which describe mixing 
of \9) and \e) due to an interaction with the environment, 
are: Ug = \s)(e\ - \0)(O\ and U s = \9)(s\ + \e)(6\. [Addi- 
tionally there is U 2 — i{\9) (e| — |e) (8\) operator belonging 
to {T 2 } — T2 representation.] These are essentially Pauli 
matrices, so the T 3 doublet may be thought of as a pseu- 
dospin (not transforming, of course, as a true spin under 
spatial rotations). 

For the Kondo coupling, one may therefore write: 



H K = UgSg + U £ S £ + U 2 S 2 



(2) 



where Sg and S £ are such combinations of the operators 
describing conduction electrons scattering from the im- 
purity ion which span the r 3 representation [for (|^) to 
be invariant]. Due to symmetry arguments, those scat- 
tering operators, spanning in general some representation 
T a , can be constructed as: 



E E ^T^k')WIN,7r)4 >/3g c k , 7r . (3) 

k,k' Pl^r 



Here the operators are defined in a standard way as 
appropriate linear combinations of operators Ck belong- 
ing to the star of the momentum vector k; k numbers 
all the stars, indexes a, f3 and 7 denote irreducible rep- 
resentations (p, q and r standing for degenerate states 
spanning these representations), (ap\\/3q, jr) are corre- 
sponding Clebsh-Gordan coefficients, and Jg (k, k') are 
the coupling constants. It should be sad that the expres- 
sion (|j|) contains, strictly speaking, infinite number of 
orbital channels. Yet only those few channels which have 
the largest coupling constants survive the scaling, the 
other ones tend to die out in the course of the renormal- 
ization group transformations J^,^,0 ■ It is this tendency 
which makes the problem tractable. The determination 
of relevant channels is a complicated problem which re- 
quires a detailed analysis of the concrete impurity struc- 
ture and which presently is far from being clear. For 
the case of the (nonmagnetic) T 3 doublet, Sg j£ is a sum 
of time-reversal partners, so that the number of orbital 



channels has to be doubled. Hence, in the simplest case, 
one arrives at the two channel model, as proposed in [ p"3"| . 
The Jahn- Teller coupling is represented by: 



Hjt = 9JT {UgQg + U e Q £ ) . 



(4) 



where {Qg,Q £ ) is a local phonon mode transforming as 
r 3 . The relevant phonons are understood to be optical 
ones and their dispersion is neglected (it is not expected 
to qualitatively modify the results). In general there is 
also a non-linear coupling (along with phonons anhar- 
monicity) . 

In order to proceed, it is instructive to formally "in- 
tegrate out" electronic degrees of freedom. This can be 
done in a standard way with the following result for the 
effective imaginary time action for the phonon mode: 



S{Q} = S^ {Q} 



(5) 



g.jT 



dg / dr 
'0 



M { e 9 \T)Qg{r)+M^\r)Q £ (T) 



where {Q} is the free action, (3 is the inverse tem- 
perature, and 

M$(r)=-(Ug, e (T)) gjT = g 

is the exact non-linear dynamic "magnetization" of the 
r 3 pseudospin. Eq. (||) contains both the electronic adia- 
batic potential for Qg t£ and a retardation term originat- 
ing from the long time relaxation of low energy electronic 
excitations. 

Consider the adiabatic potential first. It arises 
from the static magnetization (i.e. calculated for r- 
independent Qg_ £ ). As it is evident from the structure of 



and 



M, 



(9) 



(Qe,e/Q)M{gQ) with Q 2 = Q 2 + Q 2 , 
M being the transverse magnetization. Though an ex- 
plicit analytic expression for M is generally lacking, in 
limiting cases it can easily be determined. From essen- 
tially free pseudospin magnetization in high fields (tem- 
peratures) it changes to a well known logarithmic be- 
haviour in low fields (temperatures). Concerning the 
cross over, a set of numerical data for a spin isotropic 
case P3|] and analytic formulae at the Emery-Kivelson 
line H are available. The electronic part of the adiabatic 
potential adds a negative contribution to the elastic en- 
ergy (I/2)/iWg(5 2 , where /i is the mass and u>q- the bare 
frequency for the Q-mode. Hence the equilibrium posi- 
tion Q = may become unstable and shift to a finite 
Q = Qq. By differentiating the adiabatic potential in Q 
one arrives at the following equation for Qq: 



/iw 2 Qo = 9JTM(g JT Q ). 



(6) 



It is worth noticing that, since the magnetization should 
be a monotonic function of the field, this equation can 
have only one nontrivial (Qq =/= 0) solution that appears 
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at temperatures below Tjt which is determined by the 
following condition: 



HU>1 = gj T x{T JT ), 



(7) 



where x is the linear susceptibility. If the Kondo cou- 
pling (^|) is disregarded [as well as the electron-phonon 
coupling in ([!])] then a conventional Jahn- Teller effect 
takes place and Tjt is related to the energy gain due to 
the zero-temperature lattice distortion Wjt — ^jt/^M^o 
(actually Tjt = 2Wjt)- In the case when this energy is 
larger then the characteristic Kondo energy scale, Tjc, 
the Kondo coupling (even if kept finite) does not essen- 
tially influence the formation of the lattice distortion. 

One therefore arrives at the following conclusion: the 
relative importance of the Kondo- and the Jahn- Teller 
couplings is determined by the interplay of the parame- 
ters Wjt and T K ■ 

Assume that 



W JT < T K 



(8) 



(i.e. the Jahn- Teller coupling is weak). Then one can ex- 
pect Tjt <SC Tk so that the low temperature susceptibil- 
ity x(T) = (1/2ttT k ) \n(T K /T) may be used to estimate 



Tjt = T K exp {-2ir /j.u>lT K / g JT } 



(9) 



At this temperature a nonzero Qo appears first, then it 
increases with lowering the temperature reaching T = 
value which, since a magnetic field scales as y/T , is 
estimated as 



Qo 



Tr / Tjt 
gjT V T K 



(10) 



The energy gain Wjt due to this lattice distortion differs 
from the bare Wjt] it is determined by the energy scale 
(^): Wjt = Tjt/4w. It should be emphasized that for 
the Jahn- Teller effect to be fairly developed, the lattice 



distortion (10) should exceed the characteristic scale for 



zero-point fluctuations around the new minimum: 



v/((AQ) 2 > « Qo; 



(11) 



otherwise the Jahn- Teller coupling does not lead to any 
significant effects (does not really cause a distortion) and 
can be treated by perturbation theory [ p4[ . 

Theoretically, the case when both inequalities (|J) and 
( |i"T| ) are fulfilled is probably the most interesting one - 
because of its dynamics. The system can tunnel between 
different equivalent minima in the Q = Qo manifold. Be- 
yond the adiabatic approximation, each tunneling event 
causes long relaxation process involving an adjustment 
of the electronic state (the "Kondo cloud") to a new 
phonons configuration. This is somewhat similar to what 
was proposed by Yu and Anderson |2^] for ^415 com- 
pounds, though here the situation is much more compli- 
cated since a dynamic response of the two channel Kondo 



system to a rotating (or jumping by 27r/3 depending on 
the importance of warping terms) magnetic field comes 
into play ^J. One should mention that, since the lattice 
distortion arises from a competition of two quadratic in 
Q energies (elastic and electronic parts of the adiabatic 
potential) the resulting minimum is going to be rather 
flat so that the inequality ([ll]) leads to an exponential 
constraint on the mass of the Q-mode: 



V > 9jt/ t kTjt- 



(12) 



Turning now to the case of strong Jahn- Teller cou- 
pling, opposite to (H), one observes that if the lattice 
distortion is standardly given by Qo — g JT /^o then 
QjtQo is much larger then Tk. So one deals with a 
high field susceptibility when the pseudospin behaves as 
if it is free. Hence the Kondo interaction (H) is irrel- 
evant. As it is well known j2(J, the degeneracy of the 
Qo = Q\ + Qg minima manifold is lifted by phonons an- 
harmonicity and non-linear couplings (warping) . [Notice 
that the inequality ( |IT|) is still to be satisfied.] Result- 
ingly, three equivalent minima are formed at the posi- 
tions (0,1), (V3/2,l/2), (-- y/3/2,1/2) or at the posi- 
tions (1/2,^3/2), (0,-1), (-1/2.V3/2) depending on 
the details of the warping potential. With each of these 
minima a state of the system described by a product of 
the phonon and the pseudospin wave functions is asso- 
ciated. The three degenerate states will be denoted as 
|1),|2), and |3). As a matter of fact, the overall cubic 
symmetry of the total Hamiltonian ([!]) can not be vio- 
lated. This is reflected by tunneling rates for tunneling 
of the system between the states |i) (z = 1,2,3) being fi- 
nite. In the conventional theory of the Jahn- Teller effect 
p0| , p6 {27 1 , a diagonalization of the tunneling Hamilto- 
nian leads to "vibronic states": the ground doublet T3, 
spanned (in the lowest approximation in the tunneling 
rate) by the wave functions |e) = [2|1) - |2) - \3)]/VE 
and 1 9) = [\ 2) - 1 3)1 /a/2, and the excited singlet with 
the wave function |fi) = [|1) + |2) + \3)]/VH separated 
from the ground doublet by an energy AE [of the order 

of W JT exp(-3W JT /2uo)]- 

It needs to be realized that in the case of an impurity 
ion in a metal, the distorted configurations strongly inter- 
act with conduction electrons. This interaction has noth- 
ing to do with the Kondo coupling (||) - it is caused by the 
conventional electron-phonon interaction present in (Q). 
Indeed, (i) because of the latter interaction, the conduc- 
tion electrons experience different scattering potentials 
when the impurity complex is in different distorted con- 
figurations. Hence the screening effect. One may expect 
the corresponding dimensionless coupling constant to be 
of the order of unity (if the same is true for the bulk 
electron-phonon coupling). The role of the screening is 
to renormalize (decrease) the energy splitting AE. Addi- 
tionally, (ii) the electron density fluctuations modify tun- 
neling barriers between the configurations \i). This gives 
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rise to assisted tunneling processes (analogously to the 
case of two-level systems Jul). These processes, though 
in several orders of magnitude smaller then the screening 
interaction, lead to a nontrivial scaling making the prob- 
lem formally similar to the Kondo effect. As a result, one 
is dealing with a "Jahn-Teller-Kondo" phenomenon when 
the local electronic configuration, the lattice distortions 
and the conduction electrons form a coherent state. The 
coupling is of the form Ylij where Vij are the 

electron scattering operators of appropriate symmetry 
pl[ . The conduction electrons can never be disregarded. 
The present formulation gives rise to several low energy 
scales (23j . At temperatures AE <C T <C Wjt the corre- 
sponding model resembles a very anisotropic pseudospin- 
one model which critical behaviour is to be studied [^9| . 
At the lowest temperatures (T <C AE) one again expects 
a two channel behaviour (within the vibronic doublet) 
which can however not be further quenched. 

To conclude, the problem of the interplay of the 
Kondo- and the Jahn- Teller effects is formulated. Pa- 
rameters regions where the system exhibits qualitatively 
different properties are determined. The T$ ground dou- 
blet is discussed as an example though the considerations 
can also be applied, with minor adjustments pc|| , to other 
degenerate representations. I believe that various ideas 
presented here deserve further investigations. 

I am thankful to Daniel Cox for helpful discussions and 
comments. 



Gogolin, and P. Nozieres, to appear in Phys. Rev. B. 
For a discussion establishing a link between this ap- 
proach and Bet he-ansatz methods see A.M. Tsvelik, 
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